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O I Abstract. How massive were the first stars? This question is of fundamental importance for galaxy 

D . formation and cosmic reionization. Here we consider how protostellar feedback can limit the mass 

' of a forming star For this we must understand the rate at which primordial protostars accrete, how 

they and their feedback output evolve, and how this feedback interacts with the infalling matter 
We describe the accretion rate with an "isentropic accretion" model: lii^, is initially very large 
(0.03Mq yr^' when m^, = IMp) and declines as m'^^l'' . Protostellar evolution is treated with a model 
that tracks the total energy of the star. A key difference compared to previous studies is allowance 
J> ' for rotation of the infalling envelope. This leads to photospheric conditions at the star and dramatic 

CO I differences in the feedback. Two feedback mechanisms are considered: HII region breakout and 

OO ■ radiation pressure from Lyman- a and FUV photons. Radiation pressure appears to be the dominant 

' mechanism for suppressing infall, becoming dynamically important around 20 Mq . 

(N 

O THE COLLAPSE OF PRIMORDIAL GAS CLOUDS 



Recent numerical studies have followed the gravitational collapse of perturbations from 
O \ cosmological to almost stellar dimensions [[I], Baryon-dominated clouds, cooled to 

■ about 200-300 K by trace amounts of molecular hydrogen, form at the centers of dark 

. matter halos. For > 10"^ cm^-^ the cooling rate becomes independent of density, and 

so the dissipation of gravitational energy in the densest regions gradually raises the 
temperature. In the simulation of Abel et al.[[Tp, the gas cloud is centrally-concentrated 
(p o« r^^p; kp ~ 2.2) and is contracting quasi-hydrostatically with infall speeds about 
5^ I one third of the sound speed. In addition to thermal support, the cloud is filled with 

a turbulent cascade of weak shocks (T. Abel, private comm.). The structure can be 
described by an approximately hydrostatic polytrope with 7p = I + l/n = kp/kp = 
2(1 — l/fcp) = l.l, where P r^^p. The contraction is akin to the maximally sub-sonic 
Hunter^] settling solution, for which the accretion rate is a factor 0* ~ 2.6 greater than 
the classic Shu[^ solution. This accretion rate can be expressed in terms of the entropy 
parameter of the polytrope, K = P/ p'^p and the collapsed mass, M ~ m*. 



{3-kp)klK^ 



-> 0.026i^'^^/^ Mq yx-\ (1) 

V^0/ 



where j = 3(1 — 7p)/(4 — 3/^). The numerical evaluation assumes (^^ = 2.6, since jp is 
not too different from one, and = 1.88 x 10^2(^/300 K)ny'^4icgs = 1.88 x lO^^^'cgs. 



We set the temperature normalization a factor 4/3 higher than is seen in simulations [|T|, 
to allow for partial pressure support from sonic and isotropic turbulent motions; i.e. 
T is an effective temperature. The small ratio of turbulent to thermal support is in 
marked contrast to contemporary massive star formation[^. In primordial clouds it 
is the microphysics of Hj cooling that determines both the evolution (via jp) and 
normalization (via K' and 0*) of the accretion rate. Collapse assuming constant K - 
"isentropic accretion" - agrees with 1-D numerical studies[0, |]] (Fig. la). 

The mass-averaged rotation speed is a fraction, f^^^ ~ 0.5, of Keplerian, approxi- 
mately independent of radius[jT]]. Assuming angular momentum is conserved inside the 
sonic point, r^p, leads to a disk size = fl^^r.p = 3A{f^^^/0.5)^{M/MQ)'^/^K'-^^/^ A\J. 
Matter falls onto this disk at all radii r < r^^, as well as directly to the star. We follow 
Ulrich[^] in describing the density distribution of the rotating, freely-falling envelope. 

PROTOSTELLAR EVOLUTION 

At densities «^ ~ lO'^ cm an optically thick protostellar core forms[^, bounded by 
an accretion shock. The size of the protostar then changes as it accretes matter and 
radiates energy. For spherical geometry, the high accretion rates typical of primordial 
star formation lead to optically thick conditions above the accretion shock[ll^, pi]] . 
Accretion energy is advected, which swells the star. However, for collapse with angular 
momentum, a disk forms and we expect photospheric conditions over much of the stellar 
surface. We employ a polytropic stellar structure and extend the energy equation model 
of Nakano et al.[p^, that includes gravitational, ionization-dissociation, and D-burning 
energies, to allow for optical depth in front of the accretion shock. We also allow for 
some fraction of disk accretion, depending on /^ep- We model the "luminosity wave" 
expansion feature [PT)|, [TT]] as a relaxation to a more compact {n = 3) state after the star 
is older than its Kelvin time: to conserve energy the outer radius of the star grows by 
a factor, which we estimate from Omukai & Palla[pT]]. We set n = 1.5 or 3, depending 
on convective stability. During optically thick accretion, before the luminosity wave, we 
derive an effective value for n by comparing to the results of Stahler et al.[|TDp. Figure lb 
shows this case and isentropic models with f^^^ = and 0.5. The two spherical models 
undergo similar evolution, while the rotating model forms a smaller protostar once the 
photosphere retreats to the stellar surface. The protostar is then supported by D core and 
shell burning, as in present-day star formation. The large accretion rate causes the star 
to join the main sequence only at relatively high masses. The photospheric temperature 
is much hotter than in the spherical models, leading to significant fluxes of ionizing and 
FUV radiation. These form the inputs for the feedback model. 

FEEDBACK VERSUS ACCRETION 

The lifetimes of primordial massive stars converge to about 2 Myr[[T^] so eq.([T]) implies 
an upper limit to stellar masses of ~ 2000A/q. However, other feedback processes are 
likely to intervene before this. Once the flux of ionizing photons from the protostar 



(a) Evolution of Mass Accretion Rate 



(b) Protostellar and Photospheric Radii 
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(c) HII Region Geometry near Breakout 
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g (d) Mass Scales of Feedback Processes 
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FIGURE 1. (a) Protostellar accretion rate as a function of the collapsed mass (~ in these models). 
Solid line: fiducial isentropic accretion model (K' = I) from eq. (|l]); dotted[^ and dashed lines extrapo- 
lated from 1-D numerical studies; long-dashed line is m^, = 4.4 x lO^^Mg yr^i used in the protostellar 
evolution models of Stahler et al. [|o| and Omukai & Palla[pl]] . (b) Evolution of protostellar radius (lower, 
thick lines), which is the location of the accretion shock, and photospheric radius (upper, thin lines). The 
spherical constant accretion rate test case (dotted) is compared to other calculations |[T^ of (squares) 
for the evolution before H burning. Also shown are the spherical (dashed lines) and rotating (f^^ = 0.5, 
solid lines) isentropic accretion models. The initial condition is taken from 1-D hydrodynamical collapse 
simulations ||^. Note that the photospheric and protostellar radii are the same in the rotating model for 
> 0.3Mg.(c) Geometry of the Hll region (shaded) at polar breakout when r^^^^ = rg, the gravitational 
radius for the ionized gas sound speed. The protostar is at (0,0) and the disk is in the z = plane. Dashed 
streamlines show infall. (d) Mass scales of feedback processes versus the rotation parameter, f^^^- HII 
region breakout at the pole and just above the equator occur at masses traced by the lower and upper 
dashed Unes, respectively. Ly-a and FUV radiation pressure becomes greater than twice the radial infall 
ram pressure (evaluated at angle 7z/3 from the pole) for masses above the solid line. 



is greater than that of neutral H to its surface, an HII region forms. Approximating 
sectors as independent, we calculate the extent of the ionized region. Accretion may be 
suppressed if the HII region expands to distances greater than rg, where the escape speed 
equals the ionized gas sound speed, ~ lOkms^. For the fiducial model this occurs at 
the poles when = 90Mq (Fig. Ic) and at the equator when = 140Mq. In this 
calculation we assumed a free-fall density distribution. In reality the ionizing radiation 
force [[T^ decelerates and deflects the flow. For collapse with angular momentum, most 
streamlines do not come too close to the star, so the effect discussed by Omukai & 
Inutsuka[|T4|] of HII region quenching due to enhanced densities is greatly weakened. 
In fact, we anticipate that deflection is more important in reducing the concentration of 
inflowing gas near the star, so that the HII region becomes larger. 

A second feedback effect is radiation pressure from Ly-a photons created in the 
HII region and FUV photons emitted by the star. These photons are trapped by the 
Lyman series damping wings of the neutral gas infalling towards the HII region. The 
energy density builds up until the escape rate, set by diffusion in frequency as well 
as in space[0], equals the input rate. The resulting pressure acts against the infall ram 
pressure. We use the results of Neufeld[jT^ to aid our numerical calculations. Radiation 
pressure becomes greater than twice the ram pressure at ~ 20Mq for typical /j^gp. 

The enhancement of radiation pressure above the optically thin limit is by a factor ~ 10'^. 
Infall is first reversed at the poles, which would allow photons to leak out and reduce the 
pressure acting in other directions. We shall examine this scenario in a future study. 
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